Competing Magnetic Phases on a "Kagome Staircase" 
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We present thermodynamic and neutron data on NisA^Os, a spin-1 system on a kagome staircase. 
The extreme degeneracy of the kagome antiferromagnet is lifted to produce two incommensurate 
phases at finite T - one amplitude modulated, the other helical - plus a commensurate canted 
antiferromagnet for T — > 0. The H — T phase diagram is described by a model of competing first 
and second neighbor interactions with smaller anisotropic terms. Ni3V20s thus provides an elegant 
example of order from sub-leading interactions in a highly frustrated system. 

PACS numbers: 75.10.Jm, 75.25,+z, 75.30.Kz 



Geometrical magnetic frustration leads to unusual low 
temperature spin order and dynamics and presents new 
challenges for the theoretical understanding of magnetic 
systems. Frustrated materials are often characterized 
by triangle-based lattices and short-range antiferromag- 
netic (AF) interactions. 1 Of particular interest has been 
magnetism on the two-dimensional (2D) kagome lattice, 
which consists of corner-sharing triangles. While the 
Heisenberg spin- 1/2 model appears to have short range 
spin correlations and a gap to free spinons, 2,3 the S — > oo 
classical model has Neel order with a \/3 x \/3 unit cell 
at temperature T = 0. 4 Materials that approximate the 
kagome AF can be expected to lie close to a quantum 
critical point, and indeed early work on the kagome sys- 
tem SCGO exposed a spin liquid phase possessing a large 
fraction (15%) of the total spin entropy and short range 
y/S x y/3 order AS Later work on jarosite systems showed 
different "q = 0" long range order apparently favored by 
inter layer interactions. 7 

Here we study Ni 3 V 2 8 (NVO) in which the S = 
1 Ni 2+ spins form the orthorhombic kagome staircase 
structure shown in Fig. 1(a) A This structure has the co- 
ordination and two-dimensionality of the regular kagome 
lattice, but the kagome planes are buckled. The system 
is particularly attractive because its complex magnetic 
phase diagram can be understood on the basis of an em- 
bellished kagome spin hamiltonian. The model we intro- 
duce also applies to the isostructural compounds where 
Ni is replaced by Cu 9 or Co™ Although the symmetry 
of these compounds is the same as that of NVO, their 
phase diagrams are very different. As indicated below, 
this difference results from a small quantitative change in 
the parameters which dictate how frustration is relieved. 

A previous study of the magneto-thermal response in 



polycrystalline NVO revealed four zero field phase tran- 
sitions with Qw/Tn > 5, where Qw is the Weiss con- 
stant and T/v the magnetic ordering temperature™ In 
this letter we report an unexpectedly rich anisotropic 
field-temperature (H — T) phase diagram (Fig. 2), with 
high and low temperature incommensurate (IC) phases 
(HTI and LTI) and two commensurate (C and C') spin 
structures. These magnetic structures are determined via 
neutron diffraction. We also explain the salient features 
of NVO by a model, in which the spine (Ni s ) and cross-tie 
(Ni c ) spins interact via nearest neighbor (NN) and sec- 
ond nearest-neighbor (SNN) isotropic Heisenberg inter- 
actions. In addition, and consistent with crystal symme- 
try, it is necessary to take account of the Dzyaloshinskii- 
Moriya (DM) interaction and magnetic anisotropy. 

Symmetry is key to understanding the ordered phases 
that spring from the kagome critical state in NVO™ In 
the presence of AF ordering on the spine sites, isotropic 
NN interactions produce zero mean field on cross-tie sites. 
In this regard, NVO is reminiscent of Sr2Cus04Cl^ and 
of some "ladder" systems of recent interest™ However, 
the structural anisotropy of NVO induces interactions 
not usually considered in frustrated systems. First, be- 
cause the Ni06 octahedra are edge-sharing, the NN Ni-O- 
Ni bond angle is close to 90° so the NN and SNN Ni-Ni 
interactions are weak and similar in strength. Second, 
the symmetry of the crystal structure admits a DM in- 
teraction among the NN spine spins. 12 Third, anisotropic 
pseudo dipolar (PD) exchange interactions between spine 
and cross-tie spins induce both a uniform and a staggered 
moment on the cross- tie sites™ These interactions add 
to the usual isotropic NN super-exchange interaction to 
produce the observed rich H — T phase diagram. 

Single crystals of NVO were grown from a BaO-V20s 
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FIG. 1: (a) Structure of NVO, showing the cross-tie Ni c (blue 
(gray)) and spine Ni s (red (black)) sites, (b)-(c) indicate the 
spin structures in the incommensurate phases. + and - indi- 
cate spin components along b. Symbol sizes scale with the 
dipole moment, (d) indicates the symmetry of the low T com- 
mensurate spin structure. Spin canting has been exaggerated 
for clarity and the relative symbol sizes for spine and cross 
tie spins are not to scale. Subsequent layers are displaced by 
(a + b)/2 with spine spins satisfying Eq. (1). Lattice param- 
eters serve as axis length units. 

flux and powder samples were synthesized with stan- 
dard techniques i&H The uniform magnetization, M, 
was probed using a commercial SQUID magnetometer. 
The specific heat, C, was measured with a commercial 
calorimeter using the relaxation method for T > 2 K and 
the semi-adiabatic method for lower T. Powder and sin- 
gle crystal neutron diffraction measurements were car- 
ried out at the NIST Center for Neutron ResearchiAA 
The space group of NVO is Cmca (No. 64) 8 with lat- 
tice parameters a=5. 92197(3) A b=11.37213(7) A, and 
c=8. 22495(5) A at T = 1.5 K. Throughout we index 
wave vectors in the orthorhombic reciprocal lattice with 
a* = 27r/a, b* = 27r/6, and c* = 2tt/c. Representative 
specific heat data are in Fig. 3(a) for a magnetic field (H) 
of 0, 5, and 8 T along c. As in previous zero field measure- 
ments on powder samples, there are four peaks in C(T). 10 
The entropy reduction associated with these phase tran- 
sitions is determined by AS = Jq ° K (C/T)dT, after sub- 
tracting an estimate of the lattice contribution obtained 
from the non-magnetic structural analog ZnaX^Og. We 
find AS « 7.9 J/mole K, or 87% of Rln3, which is close to 
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FIG. 2: Phase diagram for NVO as a function of temperature 
and magnetic field applied along the three principal crystal- 
lographic directions. For H || c no true phase boundary sep- 
arates the P and C phases. White areas were not probed. 

that expected for ordering among spin-1 Ni 2+ ions. We 
infer that the specific heat peaks mark phase transitions 
to unique structures involving the Ni 2+ spin-1 degrees of 
freedom. The H = peaks at 2.2 K, 6.3 K, and 9.1 K 
indicate second order phase transitions, whereas the 3.9 
K peak marks a first order transition. 

Through extensive specific heat measurements, we de- 
termined the phase boundaries shown in Fig. 2. These 
were confirmed by susceptibility (x) and magnetization 
measurements (see Fig. 3(b)), which provide additional 
clues to the nature of the phases. The susceptibility ex- 
hibits significant magnetic anisotropy. As T is reduced 
and the C phase is entered, there is a sharp jump in M, up 
to 3.5% of the Ni 2+ saturation moment for H || c, which 
indicates a weak ferromagnetic (FM) moment along c. 
With H || a, there is a sharp drop in M. Finally, for 
H || b, there is no sharp feature indicating no FM mo- 
ment along b. A surprising result of this study is that 
the Tph and Thl transitions do not produce observable 
anomalies in x(T). In a field of 0.1 T the magnetization 
anomaly at Tph is less than 4 x 10 -5 /i^/Ni or 0.3% of 
the signal while it is less than 4 x 10 -6 /ie/Ni or 0.03% 
at Thl- Nonlinear susceptibility measurements likewise 
produced no indication of these phase transitions. 

Neutron diffraction, however, reveals temperature de- 
pendent magnetic Bragg peaks at Q = (2n + 1 ± 
q/a*, 2m + 1, 0) and Q=(2n + 1 ± q/a*,2m + 1, 2m + 1) 
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FIG. 3: (a) Specific heat of NVO, in zero field and for H || c. 
(b) Longitudinal magnetization versus T for H — 0.1 T 
along the three principal crystallographic directions, (c) Inte- 
grated intensity of commensurate and incommensurate mag- 
netic Bragg peaks at Q = (110) and (1 ± q/a*, 1,0), respec- 
tively, (d) Temperature dependence of the incommensurate 
magnetic wave vector. In the C and C phases we believe the 
incommensurate peak reflects a meta-stable minority phase 
as it is only present after cooling through the HTI and LTI 
phases and can be fully suppressed by field cycling. 

for Tlc < T < Tph- The peaks are resolution limited 
indicating a correlation length in excess of 500 A. The In- 
dependence of the peak intensities is shown in Fig. 3(c). 
Anomalies are apparent at the three high T transitions, 
Tph, Thl, and Tecs and the peaks vanish in the H = 
paramagnetic (P) phase. The absence of an anomaly to 
the level of 0.5% (see inset) in the T-dependence of the 
(110) magnetic Bragg peak through the phase transition 
at Tec indicates that this transition involves degrees of 
freedom that are decoupled from the prevailing AF order. 
Fig. 3(b) shows that the weak FM moment is also un- 
changed through this transition. Nonetheless we believe 
the specific heat anomaly at Tec is intrinsic as it was 
observed in all samples studied (1 powder and 5 crystals). 

The T-dependence of the characteristic magnetic wave 
vector, g, is shown in Fig. 3(d). Again there are anoma- 
lies at all the upper transitions but not at Tec- In 
phases HTI and LTI, q varies continuously indicative of 
an IC magnetic structure. The C phase is commensu- 
rate though cooling through phases HTI and LTI yields a 
metastable remnant of the IC modulation. To determine 
the spin structures in the HTI, LTI, and C phases we 
collected zero field (ZF) magnetic Bragg intensity data 
for 170 peaks in the (hkO) and (hkk) planes at T = 7 K, 
5 K, and for 70 peaks at T = 0.1 K after ZF cooling. We 



analyzed the data using group theoretical classification 
of the possible spin structures^ 

In the HTI phase, we limited consideration to magnetic 
structures that form a single irreducible representation of 
the corresponding space group, because we reject the pos- 
sibility of a multicritical point where more than one irre- 
ducible representation simultaneously become critical. 15 
Irreducible representation T^A provides an excellent ac- 
count of the HTI phase with a reliability coefficient 
R=17%. The corresponding magnetic structure is illus- 
trated in Fig. 1(b). At T = 7 K the wavelength of the a- 
modulated structure is A m = 27r/(a*— q) = 1.37(1) a with 
an amplitude vector = (1.12(4), 0.04(9), 0.01(8)) fi B 
for spine spins and = (0, —0.3(1), 0.00(6)) \i B for 
cross-tie spins. There is a phase shift of 0.4(2)7r between 
the IC waves on the two sublattices. 

The LTI phase contains an additional irreducible rep- 
resentation, T2, for the c component of spine spins 
(Fig. 1(c)). For T = 5 K the spine spin ampli- 
tudes are = (1.04(8), 0.0(1), 0.01(6)) fi B and = 
(0, -0.01(6), 0.76(7))/i B . The wavevectors for the two 
components were experimentally indistinguishable and 
the phase shift of 0.3(l)7r indicates an elliptical spiral 
in the a — c plane. The cross-tie amplitude is = 
(0,-0.8(4), 0.1(1)) ii B with a phase shift of 0.5(2)tt to 
m^. There is no detectable cross tie spin component 
associated with irreducible representation m 2 c = 
(0.0(1), 0, 0) he- The reliability coefficient was R = 22%. 

In the low temperature commensurate C phase the 
data are consistent with the spin structure shown in Fig. 
1(d) which corresponds to a mixture of representations 
Ti and Tj. While the spine sublattice is fully polarized: 
ml = (0, 0.29(8), 0) fi B and m 7 s = (2.28(6), 0, 0.1(4)) fi B , 
this is not the case for the cross tie sites where = 
(-0.26(8), 0,0) hb and m 7 c = (0, 0.3(1), 0.1(8)) fi B . In 
this fit the total magnetization along c was fixed to the 
value of 0.05 ji B per Ni atom as determined from bulk 
magnetization measurements. The reliability coefficient 
was R = 11 %. Having two active representations in C 
indicates that Tec could mark the admixture upon cool- 
ing of Ti. This scenario is, however, difficult to reconcile 
with the absence of an anomaly in the T— dependence of 
the (110) Bragg intensity (inset to Fig. 3(c)). 

We now turn to a theoretical interpretation of these 
results. First, note that the dominant AF component of 
the spine magnetization in all H = phases satisfies 

m(r) = — m(r db -b) = — m(r db -c + Sh) , (1) 

where \S\ = 0.26048(6) 8 accounts for the kagome plane 
buckling. This indicates AF interactions between neigh- 
boring spines. The spin- structure within spines is con- 
trolled by competing NN and SNN isotropic Heisenberg 
interactions denoted J\ and J2. A mean field treatment 16 
indicates that for J 2 > \Ji\/^ the spine Hamiltonian is 
minimized by a mean field spin modulation with wave 
vector, q, which satisfies cos((a* — q)a/2) = — Ji/(4J2). 
Putting aside the small T-dependence of q, we deduce 
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from the experimental value (q « 0.27a*) in the LTI and 
HTI phases that J\ ~ 2.6J2. In the presence of easy- 
axis anisotropy the highest-temperature ordered phase is 
predicted, 16 in agreement with our experiments, to be a 
longitudinally modulated phase in which the spins are 
confined to the easy a axis. If the anisotropy field Ha 
is not too large (Ha < Hi), then, as the temperature is 
lowered the longitudinally modulated phase gives way to 
one in which an additional transverse modulated compo- 
nent of spin appears, growing continuously from zero as 
the LTI phase is entered. This scenario is also consis- 
tent with our diffraction data. At still lower temperature 
the diffraction data indicate the presence of a commen- 
surate AF phase. According to mean-field theory, such 
a transition can occur for sufficiently large anisotropy, 
Ha > H2. 16 Our numerical mean field calculation^ 
show that for J1/J2 = 2.6 indeed Hi > #2, so that there 
is a range of anisotropy field for which mean field theory 
predicts the observed sequence of ZF phase transitions. 

We now discuss some of the finer details of these 
phases. From the M versus H data, extrapolated to 
H = 0, we find that in the C phase there is a weak 
ferromagnetic moment. Structural considerations show 
that the DM interaction for a single spine takes the form 

Hum = Y} D - C + (" 1 )™^ b ] ' [S(») * S(n + 1)] (2) 

n 

where n labels the spins consecutively along the spine. 
Db gives rise to a linear coupling between the staggered 
moment of the spine along a and the weak ferromagnetic 
moment of the spine along c. This weak ferromagnetic 
moment can induce a ferromagnetic moment along c on 
the cross-tie spins via isotropic Heisenberg exchange. In 
addition, such a moment on the cross-tie spins can also 
arise via PD interaction between the staggered moment 
on the spines and a uniform moment on the cross-tie 
spins. Symmetry also admits a staggered g-tensor along 
the spines, the physical origin and consequences of which 
are similar to DM interactions^ The weak ferromag- 
netism explains the absence of a phase boundary between 
the P and C phase for H| |c. In the IC phases, these inter- 



actions would give rise to modulated moments along c. 
The anisotropic interactions we invoke also generate cou- 
plings between the various IC order parameters, which 
result in weak T-dependence of the IC wave vector. 17 

Next we discuss the phase boundaries between the C 
phase and the IC phases. Barring a multicritical point, 
these must be first order transitions. For H||c the Zee- 
man energy, —HM, of the FM moment (in the C phase) 
explains why the transition temperatures Tlc an d Thc 
increase linearly with increasing H. For H _L c the 
Zeeman energy does not appear and the phase bound- 
ary of the C phase should be quadratic in H [Tn(H) = 
Tn(0) +&H 2 ] as it depends on the differences in the sus- 
ceptibilities of the phases involved. In particular, when 
H||a, the longitudinal susceptibility of the C phase is 
small and the coefficient a is negative, disfavoring the C 
phase. The other phase boundaries (Tph and Tjjl) are 
also expected to be quadratic in H and the experimental 
phase diagrams are consistent with this although for the 
HTI phase when H is along b, the coefficient a is unusu- 
ally small. This fact is linked to the absence of anomalies 
in x at the HTI phase boundaries. Both features may be 
a consequence of a frustrated and weakly connected spin 
system where phase transitions occur from a strongly cor- 
related state with short range AF order. 

In summary, we have studied the phase diagram of the 
spin-1 kagome staircase M3V2O8. We find that although 
this phase diagram is quite complicated, it can be under- 
stood on the basis of a rather simple model which reflects 
the symmetry of the crystal structure. The experiments 
and model offer a specific example of how SNN exchange, 
easy axis anisotropy, and Dzyaloshinskii-Moriya interac- 
tions can induce and control complex low temperature 
phases in a frustrated magnet. 
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